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Abstract. Various properties of single-particle Hartree-Fock ground-state solutions in TV ~ Z nuclei are 
investigated. The emphasis is on a role of single-particle mean-field in odd-even mass staggering. It is shown 
that, unlike in traditional scenario originating from the Fermi gas or macroscopic models, the symmetry energy 
contribution to odd-even mass staggering is nearly cancelled by the contribution coming from the average level 
density. It allows to construct indicators probing both pairing as well as mean-field components to the odd- 
even mass staggering. The impact of the single-particle Hartree-Fock field on Wigner energy and residual pn 
interaction in odd-odd nuclei is also discussed. 

INTRODUCTION 

Mean-field is considered as a standard nuclear model to describe medium mass and heavy nuclei. Indeed, various 
applications demonstrated ability of this relatively simple concept to describe wide range of nuclear phenomena. 
However, the accuracy of mean-field results or predictions is far from satisfactory when confronted with the present 
day, high precision experimental data. The predictions in many cases differ substantially for various parametrizations 
of the same effective interaction and the situation calls for systematic program to optimize the forces. Recent work of 
Chabanat et al. [1] can be considered as the first step towards finding comprehensive parametrization of the Skyrme 
force [2]. Indeed, the so called Lyon (SLy) forces developed in Ref. [1] are very widely used and are considered to 
be one of the best among the Skyrme forces available nowadays. 

The aim of this work is to search for certain generic features of the ground-state Skyrme-Hartree-Fock (SHF) 
solutions which may help in better understanding of basic properties of this force through the nature of self-consistent 
solutions. First section discusses a role of single-particle SHF mean-field in odd-even mass staggering (OES) which, in 
atomic nuclei, is usually attributed to pairing [3,4]. Proper understanding of mean-field and pairing contributions to 
OES is crucial to better understand and parametrize both channels. The second section investigates the possibility 
to obtain an additional binding energy in the N = Z nuclei within single-particle self-consistent SHF model or 
conventional, i.e. including only T — 1, |T Z | = 1 pairing correlations, self-consistent Hartree-Fock-Bogolyubov 
model. Finally, the third section discusses properties of self-consistent SHF solutions in time-odd channel. 



ODD-EVEN MASS STAGGERING 

The odd-even mass staggering in atomic nuclei is usually attributed to the presence of pairing correlations [3] . 
Indeed, in the simplest scenario based on standard BCS theory of superfluidity [5] , the ground state energies of three 
adjacent isotopes (isotones) can be approximately interrelated by a simple readjustment of their chemical potentials 
A: 

B(N), B{N ± 1) « B(N) + E k ± A, (1) 
where E% w A denotes energy of the lowest quasiparticle in odd-iV nucleus. Consequently, the quantity 

&i 3) (N) = ^-{B{N - I) - 2B{N) + B(N + 1)) w E k » A, (2) 
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FIGURE 1. The OES Al' as a function of neutron number N computed from theoretical masses calculated using SHF 
theory with SLy4 (upper part) and SIII (middle part) Skyrme interaction. Note that the pattern of OES is almost independent 
on the interaction. The lowest panel shows the difference of nuclear masses calculated in both cases. Unlike the OES (local 
mass correlations) calculated masses are strongly interaction dependent. 



is often interpreted as a measure of empirical pairing gap, A [6]. In Eq. (2), ir^ = (— 1)^ is the number parity and 
B(N) is the (negative) binding energy of the system with N particles. It has been noticed long time ago based 
on the Fermi gas model that the quantity (2) contains strong contribution from the symmetry energy [7]. Because 
A£, (N) is a finite step approximation of derivative of the second order: 

the symmetry energy [~ ai(A)(N — Z) 2 ] will indeed strongly contribute to OES according to this criterion. However, 
this contribution is independent on nucleon-number parity and, therefore, the symmetry energy effect can be removed 
simply by replacing the three-point indicator (2) by the higher order expression like the one containing four masses [8] : 

AW (N) ee ^ [3B(N - 1) - W(N) - B(N - 2) 

+ B(N+1)] = ^[A^(N) + A^(N-1)}. (4) 

Consequently, in the traditional picture originating from the Fermi gas or liquid drop (macroscopic) models the mean- 
field contribution to the odd-even mass staggering is predominantly due to symmetry energy and can essentially be 
removed by means of high-order indicators like (4), see [9] and refs. quoted therein. 

Contrary to the above so far commonly accepted scenario, the indicator (4) generates sizable OES when applied 
to the mass table calculated using single-particle Skyrme-Hartree-Fock model, see [4] and Fig. 1. These mass 
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FIGURE 2. The OES A^ (upper part) and A^ ; (lower part) computed from theoretical masses calculated using SHF 
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theory with SIII Skyrme interaction. Note, that rather complicated pattern of OES Aj> 4 ^ simplifies for A\f' . The A)? } shows 
clear alternating pattern with large values for even-iV and small for odd-iV. Taken from [4]. 
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calculations has been performed for nuclei with 9<Z<28 and — 2<N— Z<Q, using HFODD code (vl.75) of [10] and 
two parametrizations of the Skyrme force: SIII [11] and SLy4 [1]. The pattern of OES, which reaches from 30% 
to 50% of experimental value, appears to be generic feature of the SHF solutions i.e. is almost independent on 
paramctrization although calculated masses differ substantially for the two parametrizations employed, see Fig. 1. 

A rather complicated pattern of OES in Fig. 1 simplifies when criterion A^ 3 ) is used instead of A^ 4 \ see Fig. 2. 
Indeed, values of A^(N) are large for even-iV and small (close to zero) for odd-iV. Furthermore, it appears 
possible to explain this alternating A^ 3 ^(iV) behavior obtained directly in fully microscopic SHF calculations by 
using simple arguments based on Strutinsky energy theorem [12]. According to this theorem the results of self- 
consistent calculations can be well approximated by the microscopic-macroscopic shell-correction method i.e. the 
total binding energy can be written as B = E sp — E sp + E macl0} where 

£s P = Efe=i efe (5) 

is the shell- model energy (sum of single-particle energies of occupied states), E sp is the Strutinsky-averaged shell- 
model energy, and -E ma cro stands for the macroscopic liquid-drop energy. 
The contribution to OES from shell-model energy is: 

Se = Ag) ( A0 „ \ (1 + , N ) { e n+1 - e „) = { ^ _ ^ ![ N = ln +1 ^ 

This single-particle mechanism behind OES was first noticed in Ref. [13] and applied to OES in metal clusters [14-16]. 
The contribution coming from smooth Strutinsky energy can be expressed through average level density at the Fermi 
energy g(\): 

The average single-particle level density is g(X) = 3a/n 2 (m*/m) [13]. The empirical value of the level density 
parameter for light nuclei is afaA/8MeV what agrees well with the estimate based on realistic potentials [17]. 
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FIGURE 3. The empirical values of A^(7V = In + 1)[= A„], A\f> (N = 2n), and Ai, 4 ' (N = 2n + 1) as a function of 
neutron number. Each point represents average mean over several even-Z isotones. Note, that only (N = 2n + 1) shows 
characteristic quenching at (semi) magic gaps at N = 14, 28, 50 (no effect at N = 20 is seen) as anticipated for pairing. Thick 
gray line shows the average trend A = 12/^4 Me V. Taken from Ref. [4]. 



A correction due to effective mass is m* /m = 0.76(0.70) for SIII and SLy4 forces, respectively. Consequently, 
the contribution to OES coming from smooth Strutinsky energy equals SA^ w — 18/AMeV. Contribution from 
macroscopic energy is essentially due to symmetry energy [fj(N — Z) 2 ]. The empirical value of the symmetry energy 
strength in light nuclei is a/ = 38MeV [23]. Hence, the liquid-drop contribution to OES is SA^ « 19/AMeV what 
indeed nearly cancels out the contribution coming from the smooth Strutinsky energy. 

The alternating pattern of A^ 3 \N) allows to extract both pairing and mean-field component of OES. Indeed, at 
odd- TV mean-field component to A^ 3 ^(7V) is small and its value can be associated mainly with pairing: 



while the differences: 
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A V (N) = A<f>{N = 2n+l), 



A^ (N = 2n) - A^ 3) (N = 2n+ 1) 



(8) 



(9) 



provide information about effective Nilsson single-particle spectra [for more details concerning this aspect see dis- 
cussion and Fig. 4 in Ref. [4]]. 

Empirical values of aI. 3 -* and A^ are shown in Fig. 3. Note, that only values of A^ (N = 2n+l) show oscillatory 
behavior reflecting expected quenching of neutron pairing correlations at the magic (semi-magic) gaps at 7V=14,28, 

and 50. Rather surprisingly no quenching effect is seen at iV=20. For Al?\2n) and A^ the expected quenching 
of pairing correlations is counterbalanced by single-particle effects and no minima are seen in these cases. These 
empirical observations neatly support the above scenario. 

Similar interpretation of quantities (8) and (9) can be drown also based on seniority, equidistant-level [Richardson], 
and pairing-plus-quadrupolc models [18]. Let's consider here as an example a seniority (or pairing quasispin) model. 
The ground state energy in this model [see e.g. [6] p. 222] equals: 



B(N, s) = - — (N - s)(2Q -s-N + 2) 



where 
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where s is a seniority quantum number, stands for degeneracy of the shell and G denotes pairing strength. After 
simple algebra one gets: 



A (3) (A0 = 
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FIGURE 4. The contributions to the energy ew coming from isoscalar pairs of different angular momenta in 24 Mg (solid 
circels) and 48 Cr (open circels). The energy ew defines additional binding energy in N = Z nucleus as a displacement of its 
binding energy from the average parabolic ~ (N — Z) 2 behavior in the N / Z isostones, see inset. The details of this Ohuj 
shell-model calculation can be found in Ref. [22] . 



Note that in accordance with (8) indicator A^^N = 2n+l) probes only collective pairing energy while A( 3 )(iV = 2n) 
contains also weak, 5e — G according to (9), mean- field contribution. 

In the BCS approximation single-particle potential contributes both to A^ 3 \N = 2n) and A^ 3 \N = 2n+ 1). In 
this case one obtains Eq. (11) only in open shell regime where number of pairs is n ~ f2/2 ^> 1. In this regime, 
where BCS is considered to be very good approximation, also Abcs ~ Gfi/2[= A^(2n + l) exact }. 



WIGNER ENERGY 

It has been known since the early work of [19] that macroscopic-microscopic approaches systematically underbinds 
N w Z nuclei. Similar situation holds in semi-classical Thomas-Fermi models [20,21] as well as in spherical Skyrme- 
Hartree-Fock calculations irrespectively on parametrization of the Skyrme force [22,23]. This extra binding energy, 
which is characterized by ~ \N — Z\ type singularity at the N = Z line, is dubbed Wigner energy and is usually 
parametrized as: 

E w = W{A)\N -Z\ + d(A)ir pn 5 NZ , (12) 

where W(A) and (A) are smooth functions of mass, and Tr pn = 1(0) for odd-odd(othcr) nuclei. Parametrization (12) 
can be justified based on simple arguments of counting of np-pairs in identical spatial orbitals and on elementary 
properties of np interaction [24,9]. This estimate gives also W/d=l. In Ref. [22] (for more details see [23]) we 
constructed a family of indicators able to probe both W^(A) and d(A). It appears that indeed « d(A) 

provided that the energies of the lowest T = states are used in these indicators instead of masses of heavy N = Z 
odd-odd nuclei. It provides independent empirical argument that Wigner energy is indeed predominantly related to 
T = np-interaction. 

Unlike mean-field, the state-of-the-art shell-model calculations reproduce empirical Wigner energy very well [22]. 
One can, therefore, use this model to gain some knowledge about microscopic structure of the Wigner energy. The 
shell-model indeed relates the Wigner energy to T = interaction. It has been demonstrated first by Brenner et 
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FIGURE 5. The mass excess AM = B(;r T ) - B(x T = 1) calculated using generalized BCS plus Lipkin-Nogami theory as 
a function of T z . Different curves correspond to different values of x T being a ratio of the isoscalar to the isovector coupling 
constants. Dotted line marks the results of extended Thomas-Fermi calculations [20] while solid line without symbols Wigner 
energy according to formula from Ref. [19] . Figure taken from Ref. [26] . 



al. [25] in sd shell. The detailed microscopic analysis in sd and fp shells [22] revealed, however, rather complex 
structure of the Wigner energy in terms of nucleonic isoscalar pairs of different angular momenta. Although the 
largest contributions come from pairs coupled to J = 1 and J = J max — 5(7) for sd and fp shells, respectively, all 
matrix elements of intermediate angular momenta seem to be rather important, see Fig. 4. 

The deficiencies of conventional mean-field which allows only for T = 1, |T Z |=1 pairing can be remedied in 
generalized mean-field theory which takes into account also np-pairing correlations. It has been demonstrated in 
Ref. [26] that isoscalar np-pairing correlations can provide missing binding energy at N = Z line and its closest 
vicinity as shown in Fig. 5. The model Hamiltonian employed in [26] was very simple, based on similar pair 
counting mechanism like discussed above in connection with the Wigner energy parametrization. For such a simple 
Hamiltonian generalized BCS solutions do not mix T = and T = 1 pair correlations (see also [32]) and the energy 
gain is possible only when the isoscalar correlations are on the average stronger than the isovector. The latter cannot 
give an extra binding energy (in even-even nuclei) as long as isospin is (approximately) conserved. Simple estimates 
of isoscalar and isovector pair correlations in simple nuclei like 42 Sc suggest indeed that isoscalar correlations are on 
the average stronger than isovector, see discussion in [27,28]. However, connection between empirical interactions 
deduced from the spectra of simple nuclei and the mean-field effective pairing interaction is not obvious. Fast 
disappearance of isoparing with increasing T z = (N — Z)/2 shows up naturally in generalized BCS model and can 
be easily understood in terms of blocking of np-pairing due to proton or neutron excess [26,29]. One should mention 
here, that recent generalized BCS calculations with G-matrix interaction [31] in fpg shell (^4 ^80) do allow for 
mixing of T = and T = 1 pair correlations similar to the shell-model solutions [32]. 

The np-pairing is not necessarily the only contribution to the Wigner energy. In principle also self-consistent mean- 
field itself can contribute to the extra binding in TV = Z nuclei because of congruent nodal structure of a neutron 
and proton wave functions which can lead to stronger np-interaction as suggested by [21]. The congruence energy 
can manifest itself only in fully microscopic models. Spherical, conventional Hartrcc-Fock-Bogolyubov calculations 
with T = l,\T z \ = 1 pairing only do not show any congruence effects as shown in Fig. 6, see also Refs. [22,23]. 
This result neither depends on Skyrme interaction nor pairing force. However, deformed SHF-plus-BCS calculations 
do show around ^30% of the empirical Wigner energy. Both spherical and deformed calculations are illustrated in 
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FIGURE 6. The values of 8V pn (see Eq. (13)) calculated using theoretical masses computed by means of spherical SHF 
model with SLy4 Skyrme force and density dependent, surface active delta interaction in pairing channel (upper part) and 
theoretical masses computed by Tajima et al. [30] using deformed SHF-plus-BCS theory with SIII Skyrme interaction (lower 
part). Figure taken from [23]. 



Fig. 6 which shows values of the so called double-difference indicator [35] 

SV np = [B(N, Z) - B(N - 2, Z) - B(N, Z - 2) + B(N - 2, Z - 2)]/4 » 8 gf^'f* ■ (13) 

This quantity is particularly convenient to probe Wigner energy [22]. In deformed calculations (only even-even 
nuclei) masses computed by Tajima et al. [30] were used to calculate SV np . Very similar effect (about 30% of 
experimental value) was obtained also in our single-particle deformed SHF calculations [23]. In these calculations 
also odd-odd nuclei has been included. Rather surprisingly, an enhancement of 8V np in odd-odd N = Z nuclei comes 
entirely from the time-odd isoscalar field, see next section for greater detail. 



THE TIME-ODD MEAN-FIELDS 

Investigation of odd and odd-odd nuclei require the time-odd mean-fields to be systematically taken into account 
in the calculations. The time-odd part of the energy density T-[ odd reads: 

H° dd = C s t s 2 t + Cf s s t ■ As t + Cfs t ■ T t + C{j 2 t C7 j s t ■ (V X jt) (14) 

where s,j, T are time-odd spin, momentum, and vector kinetic densities, respectively, and subscript i[=0,l] denotes 
isospin [33,10]. 

The effective contribution of the time-odd mean-fields to the total energy E to can be investigated by comparing 
full SHF calculations with the calculations restricted to time-even fields only. The result is shown in Fig. 7. As seen 
from the figure everywhere beyond N = Z line E to is small but force dependent. Indeed, the time-odd fields are 
effectively attractive for SLy4 but repulsive for SIII forces, respectively. Also beyond N = Z line simple additivity 
for the averages [see inset in Fig. 7]: 

E to (N - Z = 2n) 7aE to {N - Z = 2n-l) + E to (N - Z = 2n + l), (15) 
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FIGURE 7. The effective energeticai contribution to nuciear mass due to time-odd mean-fieid as a function of N — Z. The 
caicuiations were performed using SHF model with Sly4 (upper part) and SIII (lower part) Skyrme forces, respectively. An 
inset represents configuration independent contributions i.e. arithmetic averages over T z =const nuclei. 



is rather well fulfilled, reflecting single-particle nature of the effect. The most remarkable observation is an enhance- 
ment of the time-odd effects in N = Z nuclei. The effect is essentially due to cancellation of two large components in 
(14), namely a repulsive CqSq field is overbalanced by an attractive C As So • Aso field. The effect reflects most likely 
spontaneous spin polarization at the nuclear surface but its nature and consequences are not yet fully recognized 
and are under study [36]. It is noteworthy, that any observable consequences of this spin collectivity may help to 
establish certain constraints on time-odd coupling constants Cq and C^ s . Note that C t s and C t As are the only 
free time-odd coupling constants for Skyrme force. The remaining coupling constants (14) are related to time-even 
coupling constants because of local gauge invariance of the Skyrme force [34] . 

It is interesting to observe also strong instabilities in time-odd fields for N ^ Z 43 Sc, 44 Sc nuclei and in their 
isobaric analogs 43 Ti, 44 V. These instabilities show up only for SLy4 force and only at the very bottom of the fan 
shell. This is, most likely, due to the extended spatial dimensions of these nuclei and due to strong attractiveness of 
Cq component of the SLy4 force at low densities where Cq(p — ► 0) = — 207.8 MeVfm 3 . For most of the commonly 
used Skyrme forces C${p — > 0) > 0, see table I in Ref. [33]. 

Let's finally consider briefly an indicator 

e pn = (-l)( A+1 )[-B(iV - 1, Z- 1) + 2B(N, Z - 1) - B(N + 1,Z - 1) + 
2B(N -1,Z)- 4B(N, Z) + 2B(N + 1, Z) 
—B(N - 1, Z + 1) + 2B(N, Z+l)- B(N + 1,Z+ l)]/4, (16) 

which probes residual pn-interaction in odd-odd nuclei, see [9]. The data indicate, that for light N ~ Z nuclei e pn 
is weekly dependent on nuclear mass and e pn ~ 500 keV. Only in the closest vicinity of N = Z line e pn is enhanced 
due to the Wigner energy, see [22,23]. Applying the indicator (16) to single-particle SHF mass table gives no effect 
i.e. e pn ~ in odd-odd N ^ Z nuclei. The cancellation of time-odd effects in (16) is a consequence of the above 



mentioned, see Eq. (15), energetical additivity of time-odd contributions in odd and odd-odd nuclei. Only in N = Z 
odd-odd cases e pn ^ as a result of an enhancement of the time-odd fields in these nuclei. It seems therefore justified 
to conclude that the residual neutron-proton interaction in odd-odd nuclei as seen through the indicator (16) goes 
entirely beyond the mean-field. 

SUMMARY 

Selected properties of single-particle self-consistent Skyrmc mean-field have been analyzed. It has been shown by 
direct microscopic calculations that the mean-field component to OES according to the criterion (2) shows strong 

nucleon number parity ttn(Z) dependence with AV being large (small) for ttn(z) = +1(-1)> respectively. This 
alternating pattern is due to (i) cancellation of the contributions from symmetry energy and average level density, 
and (ii) due to nuclear Jahn- Teller effect. It allows to interpret A^^(iV = 2n + 1) directly as due to (almost) pure 

pairing and construct higher order indicator (9) to probe also the mean-field component to OES. This scenario is 
entirely different than the commonly accepted scenario derived from the Fermi gas model where mean-field component 
to OES is due to symmetry energy and thus is independent on the number parity. 

It has been demonstrated that spherical HFB calculations including only T = 1,|T Z = 1 pairing generate no 
congruence energy i.e. VK(A) w 0. It appears also, that deformed single-particle SHF (or SHF-plus-BCS) mean-field 
can generate only relatively small fraction ^30% of the empirical Wigner energy strength VL(A) in even-even N = Z 
nuclei. Therefore, most of the Wigner energy strength in N = Z nuclei seems to be beyond the Skyrme mean-field 
and, most likely, is due to neutron-proton pairing. 

One of the most spectacular observations is strong enhancement of the time-odd effects in N — Z nuclei. Closer 
examination shows that this effect depends essentially on the first two components (more precisely is due to their 
isoscalar parts) in (14). Because Cq and C^ s are independent (free) time-odd coupling constants for Skyrme force 
any observable consequences of this effect can help to establish empirical constraints for their values. In contrast to 
the N = Z nuclei, in N ^ Z nuclei effective contributions to the total binding energy coming from the time-odd 
fields reflect simple additivity pattern between odd and odd-odd nuclei typical for a single-particle model. As a 
consequence, binding energy between valence proton and neutron in odd-odd nuclei meassured by means of the 
indicator (16) is zero. In other words e p „ is entirely beyond SHF field. 
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